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The two-photon ls 2 2s2p 3 Po —> Is 2 2s 2 x So transition in berylliumlike ions is theoretically inves¬ 
tigated within a full relativistic framework and a second-order perturbation theory. We focus our 
analysis on how electron correlation, as well as the negative-energy spectrum can affect the forbidden 
E1M1 decay rate. For this purpose we include the electronic correlation by an effective potential 
and within an active-electron model. Due to its experimental interest, evaluation of decay rates are 
performed for berylliumlike xenon and uranium. We find that the negative-energy contribution can 
be neglected in the present decay rate. On the other hand, if contributions of electronic correlation 
are not carefully taken into account, it may change the lifetime of the metastable state by 20%. By 
performing a full-relativistic jj-coupling calculation, we found discrepancies for the decay rate of an 
order of 2 compared to non-relativistic LS'-coupling calculations, for the selected heavy ions. 

PACS numbers: 32.80.Wr, 31.30.jd, 32.70.Cs 


I. INTRODUCTION 

Two-photon decay has been studied several times since 
it was originally discussed by Goppert-Mayer pQ. In low- 
Z atomic systems, the 2s — Is transitions in hydrogen¬ 
like and heliumlike ions occur primarily by two electric 
dipole photons (E1E1), and the respective decay rates 
provided by theory and experiment are in good agree¬ 
ment. These works focused not only on the total and en¬ 
ergy differential decay rates usi, but also on the angular 
and polarization correlations of the emitted two photons 
[5H9]. Detailed analysis of these two-photon properties 
have been used to reveal unique information about elec¬ 
tron densities in astrophysical plasmas and thermal x-ray 
sources, as well as highly precise values of physical con¬ 
stants m • The study of two-photon decay in high-Z 
ions also provided a sensitive tool for exploring the rela¬ 
tivistic and quantum electrodynamic (QED) effects that 
occurs in the strong atomic fields of those systems. As 
in the case of low-Z ions, predictions for two-photon de¬ 
cay rates are in good agreement with experimental data 

hthhi. 

Scarce investigations have been performed so far for 
other atomic systems with more than two electrons. In 
the case of lithiumlike ions, this lack of research might 
be attributed to almost all two-photon transitions be¬ 
ing in direct competition with dominant allowed (single 
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El) transitions, thus reducing the importance of the for¬ 
mer process in practical applications. However, this is 
not the case for berylliumlike ions with zero nuclear spin 
(/ = 0). Owing to the 0 —)> 0 selection rule, the first 
excited state ls 2 2s2p 3 Pq is metastable and its transi¬ 
tion to the ground state Is 2 2s 2 1 So is strictly forbidden 
for all single-photon multipole modes. The most domi¬ 
nant decay process is a rare two-photon transition with 
a magnetic dipole mode (ElMl) that is very sensitive to 
relativistic and electronic correlation effects and can have 
lifetimes from few decades to few minutes, depending on 
the atomic electromagnetic field of the nucleus. 

Knowledge of metastable decay rates are essential in 
collision-radiative modelling of astrophysical low-density 
plasmas that occurs in stellar coronae HE thus many 
studies have been dedicated to the measurement and cal¬ 
culation of higher-order (Ml, E 2) and hyperfine-induced 
El transitions modes [T6MT8] . First measurements of the 
metastable hyperfine-induced decay rate in N 3+ was first 
performed at the Hubble Space Telescope with important 
implications to the isotopic abundance in an observed 
nebula m ■ Values of El Ml decay rate in berylliumlike 
sulphur can also play an important role, specially be¬ 
cause the majority of stable isotopes ( 32 S and 34 S) con¬ 
tains 7 = 0 and have observable quantities in the solar 
coronae [T8l f20l f2Tj. 

Besides this astrophysical interest, there is also moti¬ 
vation for calculating the E1M1 two-photon decay mode 
coming from experiments aimed to test the standard 
model via the observation of parity nonconservation in 
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FIG. 1. (Color online) Energy atomic structure of berylli¬ 
umlike uranium relative to the final state 2s 2 x So. 


berylliumlike uranium [22] [23]. Moreover, some subtle 
X-ray lines coming from an electron cyclotron resonance 
(ECR) plasma might be attribute to charge state mecha¬ 
nisms involving the Be-like metastable state ls 2 2s2p 3 Po 
[24]. 

There are no experimental results for the E1M1 de¬ 
cay rates in zero-spin berylliumlike ions available. Only 
recently, several dielectronic recombination resonances 
were clearly identified as coming from a parent ls 2 2s2p 
3 P 0 metastable state in xenon ( 136 Xe +5 ° with 7 = 0), 
which will lead to a forthcoming measurement of the re¬ 
spective E1M1 decay rate [25] [26]. In these recently 
published works, the need of full relativistic calculations 
for this decay rate is emphasised. 

In isotopes with non-zero I the importance of the 
E1M1 mode is reduced: the hyperfine mixing between 
the term 3 Pq and the closely-lying above term 3 Pi pro¬ 
duces states with total angular momentum F / 0, thus 
circumventing the 0 —)> 0 selection rule. This drastically 
reduces the lifetime of the metastable ls 2 2s2p 3 Po state 
since it opens an FT single-photon channel. Decay rates 
for this hyperfine-induced El mode have been known the¬ 
oretically for some years HUE3ES1 and the first mea¬ 
surements have already been performed recently, both in 
laboratory [Sol EI] and in galaxy nebula m . A review 
of this topic can be found in Ref. [32 ^. 

From the theoretical point of view, the calculation of 
this two-photon decay rate offers a challenge not only be¬ 
cause berylliumlike ions have a compact electron struc¬ 
ture, which makes electronic correlation of paramount 
importance, but also due to relativistic effects, such 
as the negative-energies. Previous studies about two- 
photon decay with a Ml component have shown that this 
negative-energy contribution is mandatory for both low- 
Z and high-Z ions [33, 34] and improves gauge invariance 
[35H37] . Furthermore, similar investigations have con¬ 
cluded that the inclusion of negative-energy continuum 


gives better agreement with experimental data |38l . Both 
effects have to be efficiently incorporated in the second- 
order summation over the intermediate states that char¬ 
acterizes two-photon transitions. 

Figure [l] illustrates the compact atomic structure in 
berylliumlike uranium, where the initial, intermediate 
and final states are plotted. 

Up to now, only two estimations of the E1M1 de¬ 
cay rate for berylliumlike systems are available [39, 140] . 
both assuming a non-relativistic approximation and us¬ 
ing LS'-coupling, which for high-Z ions may lead to sig¬ 
nificant deviations. Moreover, the summation over the 
intermediate states was only restricted to the first terms, 
ls 2 2s2p 3 Pi and ls 2 2s2p ^-Pi. 

In this work, we calculate the two-photon decay rate 
of the metastable ls 2 2s2p 3 Pq state in berylliumlike ions 
considering a relativistic evaluation of the second-order 
summation in a jj-coupling active-electron. Negative en¬ 
ergies are thus included and investigated. In order to take 
into account the electronic correlation, we perform the 
evaluation of the second-order summation via a finite- 
basis-set and an effective local potential, with a few key 
intermediate states calculated using the MultiConfigura- 
tion Dirac-Fock (MCDF) method. For these evaluations, 
we consider xenon and uranium, following the reason¬ 
ing above. For elements below xenon, we notice that the 
strong electronic correlation prevents the present method 
of retrieving a reliable decay rate. A model beyond the 
active electron model is currently under investigation. 

II. THEORY 

The evaluation of two-photon related quantities have 
been discussed several times in the literature mmm, 
we, therefore, present here only the final form suitable 
for further discussion of the influence of the relativistic 
and electronic correlation effects. 

Two-photon processes are evaluated following a 
second-order perturbation theory, which overall contains 
a summation over the complete spectrum of a given 
Hamiltonian. Its elements are often referred as inter¬ 
mediate states. For the present case of the E1M1 
two-photon decay between the states ls 2 2s2p 3 Po and 
Is 2 2s 2 (terms are given for state identification), the 

differential decay rate is given by (atomic units), 


£ [S jn ( 2 , 1 ) +Si” ( 1 , 2 )] , 

jn = 1/2 

(1) 

where uq and Coq are the energies of the two emitted pho¬ 
tons, c is the light speed and j n is the total angular 
momentum of the active electron performing the tran¬ 
sition. From now on, we write configurations without 
Is 2 for shortness. The sum of both photon energies is 


dW UJ 1 UJ 2 

duoi (27r) 3 c 2 
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equal to transition energy due to energy conservation, 
E 2s2p 3 P 0 - E 2 S 2 i So = + ^2 = Wt- The two-photon 

amplitudes S' Jn (2,l) and S^ n ( 1,2) contains the summa¬ 
tion over the reduced matrix elements of the El and Ml 
multipole components, which are given by, 

5^ (2,1) = 

(2S 1 So | |-R2 | | n vjv) ( n vjv | l-Rl | | 2^ 3 P 0 ) 

„ E n~ E 2 s2 p3P 0 + W 1 

with the multipole components, electric dipole and mag¬ 
netic dipole being given by the relativistic radiative oper¬ 
ators R\ = E\ and R 2 = Mi, respectively m- sh 1 , 2 ) 
is given by an equation similar to Eq. © by interchanging 
1 with 2. In Fig. [l] are represented the first states of the 
summation for the four two-photon amplitudes allowed 
by selection rules, which are 2 < s3s 1 / 2 1 So, 2s3d 3 / 2 1 Di, 
2s2pi/2 3 Pi and 2s2p 3/2 1 P 1 , for 5 1/2 (2,1), 5' 3 / 2 (2,1), 
5' 1 / 2 (1,2) and *S' 3 / 2 (1,2), respectively. 

In this work, we consider the active electron model 
(AEM) [41], i.e, only intermediate states with variations 
of the active electron’s quantum numbers n'l'-, that par¬ 
ticipates in the transition 2 p —>> n'l'-, —>• 2s are taken into 
account in the summation over the intermediate states. 
Other intermediate states with excitation of the specta¬ 
tor electron, (like the Is and 2s occupied orbitals) are 
thus not taken into account. Atomic states are usually 
given as a linear combination of configurations within a 
MCDF or configuration interaction (Cl). We hereby de¬ 
fine a state with major contribution of a configuration 
with a spectator-orbital excitation as C exc . (7 non " exc are 
usual states within the AEM, where the major contri¬ 
bution addresses to non-excitation configurations of the 
spectator-orbital. In order to better justify the AEM, 
we give in Figure [2] a pictorial representation of one C exc 
and one (7 non - exc . The path (a) corresponds to a state 
Cnon-exc - m ^he AEM. This path is also represented in 
Fig- HI Because the radiative operator is a one-body op¬ 
erator, (7 exc states give non-null matrix elements only by 
considering either of these two cases: 

• Path (b)- Electron orbitals are almost orthogonal 
between all states, thus if the radiative operator 
connects the active orbitals, there is a small con¬ 
tribution of C exc due to (2s im |3s mt ) / 0 and 
(3s mt |2s fin ) / 0, where 2s im , 3s mt and 2s fin are 
spectator orbitals in the initial, C exc and final 
states. 

• Path (c)- The final and initial states can have a 
reasonable contribution of a configuration with the 
same spectator orbital as C exc due to configuration 
mixing. The configuration coefficients can be ob¬ 
tained either by MCDF or Cl. 

These two cases show how multiconfiguration and fully 
relaxed orbitals can play a role in two-photon processes 
by allowing C exc states beyond the AEM. 



FIG. 2. (Color online) Representation of C non ' exc and C exc 
possible intermediate states. The path (a) (solid arrows) cor¬ 
responds to a intermediate state (7 non-exc within the AEM. 
The path (b) (dashed arrows) connects an intermediate state 
C exc to the final and initial states via the nonorthogonality of 
the spectator electrons. Path (c) (dot arrows) connects cor¬ 
relation configurations of the initial and final states to C exc . 
The dashed line means a link between the configurations with 
(co) higher and (ci) minor contributions in the jj-expansion. 


For the present elements of xenon and uranium, MCDF 
calculations addressed to the initial, intermediate and fi¬ 
nal states with fully relaxed orbitals shows data that jus¬ 
tifies the use of the AEM: First, due to the strong field of 
the nucleus, even for xenon, the obtained radial orbitals 
of the spectator electrons are reasonably orthogonal in 
all initial, final and intermediate states (residues of 2%); 
Second, the jj —expansion of both initial, intermediate 
and final states are well represented by a single config¬ 
uration (all other configuration coefficients adds up to 
2 %). 

For low-Z ions or neutral beryllium, on the other hand, 
strong electron correlation does not allow the application 
of AEM. The intermediate state summation can be done 
either via an inhomogeneous four-electron Dirac Hamil¬ 
tonian (as in Ref. [12] for heliumlike ions), or by the 
introduction of C exc , which requires a careful analysis of 
the configuration mixing coefficients and orthogonality. 

In the present AEM, the evaluation of the two-photon 
amplitudes is performed by applying the finite-basis-set 
(FBS) method to the representation of the C non " exc inter¬ 
mediate states, which are eigenstates of the Dirac many- 
electron Hamiltonian. A B-spline basis set @21 mg is 
considered for a cavity of radius 60 atomic units and 50 
positive-energy and 50 negative energy states. Since the 
AEM is employed here, the FBS spectra addresses the 
active electron. A degree of correlation is introduced in 
order to match the respective orbitals obtained by the 
MCDF method. The local electrostatic potential formed 
by the Is and 2s spectator orbitals, 2uo(ls,r) +fo(2s,r) 
0, is considered in all active states, where vq (y,r) is 
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TABLE I. Optimal values of k\ and energy differences (eV) 
obtained by the FBS method without k\ optimization (Fp BS ), 
as well as the respective ones obtained by the MCDF method 
(Fmcdf). Energy of 2s2p 3 Po is relative to final state 2s 2 
x So (^t), while the rest are relative to the initial state, E n — 
^2s2 Pl/2 3 p 0 - Values provided by Ref. |l6] are also listed. 


Xe 50+ 3 P 0 - 1 So 

ki 

0.79 

7TI* 

-C'FBS 

118.5 

Fmcdf 

104.1 

Ref. [46] 
104.5 

3 Pi- 3 Po 

0.64 

0.0 

23.9 

22.8 

^i-^o 

0.62 

400.8 

430.1 

428.3 


U 88 + 

3 Po- 1 S 0 

0.69 

252.7 

258.1 

258.3 


3 Pi- 3 Po 

0.58 

0.0 

41.6 

39.9 


‘Pi-^o 

0.72 

4259.3 

4245.3 

4243.3 


given by 

vo(v,r) = J (Pv{r') 2 +Q v {r') 2 ) 2-d r ' , (3) 


Here, P v and Q v are the large and small components 
of the radial wavefunctions of a spectator orbital v and 
r> = max(r, r'). A comparison of the spectator orbitals 
of all states obtained by MCDF shows 5% differences, 
which for the evaluation of the electrostatic potential can 
be neglected. The spectator orbitals of the 2s2p 3 Po state 
are chosen for vo(ls,r) and vo(2s,r). A local statistical- 
exchange potential is also included in order to approxi¬ 
mate the non-local part of the Dirac-Fock equation. We 
follow the original procedure of Cowan [45] that defines 
this local potential for an orbital v as, 


fexc {v,r) = 


x 



Ar) 

p'{r) + 0.5/(n„ 
( 24p(r)\ 1/3 
\ * ) 


K)\ 


(4) 


where p is the many-electron total electron density and 
p'(r) is the modified total density without the contribu¬ 
tion of the v orbital, i.e., p'(r) = p(r) — min(2, e v )p u (r), 
with p u (r) being the electron density of the orbital v. 
The quantity e v is the number of equivalent electrons at 
the orbital v with principal quantum number and orbital 
angular mometum, n v and C, respectively. The func¬ 
tion (f>(r) takes into account the different influence of the 
centrifugal potential to the various orbitals as described 
in Ref. [45]. All the present wavefunctions and densities 
necessary for calculating Eqs. © and Q were obtained 
by the MCDF method. 

Next, we identify the intermediate states with the most 
relevant weight to the summations and calculate their 
most accurate MCDF energies E n . These intermediate 
states are depicted in Fig. [l] While the parameter k\ is 
set to 0.7 in Ref. [45] as the best empirical guess for the 
exchange potential, we here consider it as a free param¬ 
eter. Optimal values of this parameter are obtained by 
comparing the values of the transition energy (uj t ) and 


the energy differences E n — E 2s 2 Pl/2 3 p 0 (denominators 
of Eq. ©), obtained by the FBS method and with the 
respective ones of the MCDF method. Table [i] lists the 
optimal values of k\ that minimizes the differences be¬ 
tween the FBS and MCDF of the mentioned energy dif¬ 
ferences. The MCDF calculations were performed using 
the general relativistic MCDF code (MDFGME) [47l . 

Calculations of the decay rate were performed in both 
length and velocity gauges. The quality of the evaluation 
of the two-photon amplitudes, if the potential remains 
local in all states, is directly connected to the gauge in¬ 
variance HUES]. Although we introduced different local- 
exchange potentials in the states and MCDF energies, we 
notice that the gauge invariance is still at a level of few 
percent. 

With the application of the present formalism to the 
decay of ls2p 3 Pq to the ground state in heliumlike ions, 
and with an effective potential of fo(l<s, r), we reproduce 
the results of Ref. [48] within the respective accuracy. 


III. RESULTS AND DISCUSSION 


The results of our calculations for the 2s2p 3 Po -A 
2s 2 1 So E1M1 decay rate W are presented in Table |n| 
The obtained lifetimes corresponds to ^3 min and 12 s for 
Xe 50+ and U 88+ ions, respectively. Other allowed higher- 
order multipole contributions to this transition, like the 
E2M2 or P3M3, are severely reduced. The obtained 
value for the E2M2 decay rate in beryliumlike uranium 
of 8.2 xlO -18 s -1 shows the minimal impact to the total 
decay rate. 

Calculations were performed in both velocity and 
length gauges, showing differences of up to 4% due to the 


TABLE II. Decay rate (s _1 ) for 2s2p 3 P 0 -A 2s 2 % E\M\ 
transition in xenon and uranium. Relativistic calculations 
have been performed in velocity (V) and length (L) gauges 
for several cases: with (W opt ) and without k \-optimization 
(W non " opt ); with the summation carried without negative en¬ 
ergies (1V + ); having the energies provided by Ref. [46] (IT*); 
without the effective exchange potential (R/ no " exc ). Values of 
Refs [40) and 25] are listed. 



w° pt 

W° pt 

j^non-opt 

^non-opt 

Xe 50+ 

4.78 x 10“ 3 

4.97 x 10“ 3 4.05 x 10“ 3 

4.05 x 10 -3 

U 88 + 

8.04 x 10“ 2 

8.06 x 10“ 2 9.35 x 10“ 2 

9.35 x 10“ 2 


W+ 

W+ 


W£ 

Xe 50+ 

4.78 x 10“ 3 

4.98 x 10“ 3 5.20 x 10“ 3 

5.40 x 10 -3 

U 88 + 

8.08 x 10“ 2 

8.11 x 10“ 2 8.18 x 10“ 2 

8.20 x 10“ 2 


j^no-exc 

Ref. 20 

Ref. [2Sfl 


Xe 50+ 

5.30 x 10“ 3 

3.4 x 10“ 2 

5.2 x ltT* 


U 88 + 

8.31 x 10“ 2 

2.6 x 10 1 

4.9 x 10 1 



a Extension of Ref. m having the energy splitting 3 Po — 3 Pi 
into account. 
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FIG. 3. (Color online) Normalized differential decay rate 
for Xe 50+ and U 88+ . Values in a LS -coupling scheme are 
from Ref. [39]. The quantity y represented the energy sharing 
between both photons, i.e., y — uoi/uot . Fx figure 


different local-exchange potentials in the states. The case 
without these effective exchange potentials (VF no-exc ) re¬ 
sults in a gauge invariance of 10“ 10 %. 

Differences between the values of the decay rate with 
and without k \-optimization in Table [IT] are mostly due 
to the respective transition energies, for which the decay 
rate depends quadratically, as well on the different 3 Pi 
and 1 Pi energies. These values can also can be compared 
with the case of not considering the effective exchange 
potential of Eq. 0- Differences of up to 20% and 16% in 
xenon and uranium, respectively, shows how sensitive the 
decay of this transition is to the electronic correlation, in 
particular to the non-local part of the electron-electron 
interaction. 

Residual differences of 0.5 — 2 eV between MCDF en¬ 
ergy values and those of Ref. [46 results in relative dif¬ 
ferences of 8% in the decay rate. Most of the experi¬ 
mental observations pa S3 ED] and theoretical calcula¬ 
tions [46l f5T1 [54] of these energies are included in a energy 
range of 2 eV, resulting in differences up to 10%. 

To be conservative, we consider the uncertainty in the 
decay rate as the combined uncertainty of the previous 
effects. The final result of the E1M1 decay rate is thus 
equal to (5 ± 1) x 10 -3 s -1 and (8 =b 1) x 10 -2 s -1 for 
xenon and uranium, respectively. 

In contrast to previous studies of the negative con¬ 
tinuum, where it shows that its contribution has to be 
included in relativistic calculations of two-photon decay 
rates [33] [34], the present case is of order of few per¬ 
cent, even for berylliumlike uranium ions. Following the 
semirelativistic approach of Ref. [48] . the estimation of 
the negative-continuum contribution to the decay is pro¬ 
portional to cj 3 /Z 2 . Previous studies deal with transi¬ 
tions between principal quantum numbers (e.g. [34]). 

for which the transitions energies scales as Z 4 . In the 
present case, the transition addresses the same quantum 
number and scales roughly as Z. This makes a smaller 
contribution of the negative energy continuum. 

We notice evident differences relative to previous cal¬ 
culations by factors from 10 to 300. The differences can 
be attributed to our full relativistic approach in a jj- 
coupling scheme. This can be further investigated in the 


differential decay rate that is illustrated in Fig. [3] where 
it is shown normalized values (to the integral) of Eq. 0 - 
Here, values for Xe 50+ and U 88+ ions obtained in this 
work and by Ref. [39] are displayed, which show evident 
differences in the differential decay rate. 

The values of Refs. [39] 00] where obtained by con¬ 
sidering only the 2s2p 3 Pi and 2s2p x Pi states in the 
intermediate-state summation and were calculated in a 
non-relativistic PP-coupling framework. Moreover, the 
non-relativistic form of the electric and magnetic dipole 
operators was also employed in Refs. [39], 40], which 
forbids intercombination transitions with a spin-flip of 
the total spin in a PP-coupling. Therefore, spin-orbit 
and spin-spin interactions were included in first approx¬ 
imation in order to mix the 3 Pi and l P\ terms. For 
highly charged ions, intercombination transitions are al¬ 
lowed in a j j-coupling scheme with relativistic wavefunc- 
tions, as the spin-orbit interaction is already included 
non-perturbately. Other investigations of the P1P1 have 
shown that relativistic effects increase the decay rate by 
30% H3ES] in heliumlike Xe. In the present case, the Ml 
mode is even more sensitive to the PP-coupling scheme 
that is not appropriate for highly charged ions, where the 
strong spin-orbit interaction is included perturbately. A 
similar factor of 300 was already obtained in a relativistic 
calculation [56] . 

IV. CONCLUSION 

We have presented the results of the two-photon for¬ 
bidden El Ml decay rate for two selected heavy elements 
obtained with an effective potential. The limitations of 
the active electron model for this particular decay is in¬ 
vestigated and found that while this approach cannot 
be applied to low- and middle-Z ions, for berylliumlike 
Xe and heavier elements, each state is well described by 
a single configuration with orthogonal orbitals. There¬ 
fore, excitations of the spectator electron that forbids 
the use of this model can be neglected. We have found 
a negligible contribution of negative-energy states to this 
decay rate, which is in agreement with semirelativistic 
estimations. On the other hand, we observe significant 
relativistic effects relative to non-relativistic calculations 
performed for middle -Z ions, which can be attributed to 
the fact that the LP-coupling scheme is not appropriate 
for the evaluation of this decay rate in highly charged 
ions. 
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